Finite-range effects in dilute Fermi gases at unitarity 
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We develop a theoretical method going beyond the contact-interaction approximation frequently 
used in mean-field theories of many-fermion systems, based on the low-energy T-matrix of the pair 
potential to rigorously define the effective radius of the interaction. One of the main consequences 
of our approach is the possibility to investigate finite-density effects, which are outside the range 
of validity of approximations based on 5-like potentials. We apply our method to the calculation 
of density dependent properties of an ultracold gas of ''Li atoms at unitarity, whose two-body 
interaction potential is calculated using ab initio quantum chemistry methods. We find that density 
effects will be significant in ultracold gases with densities one order of magnitude higher than those 
attained in current experiments. 



I. INTRODUCTION 

In recent years, remarkable experimental advances in 
the field of trapped alkali gases made it possible to tune 
the interaction of fermionic atoms and study experi- 
mentally the transition betvifeen a BCS state of paired 
fermions to a BEG state of diatomic molecules [1]. 

Several interesting phenomena occurring in these and 
other fermionic systems are usually analyzed by using 
a mean-field treatment of a model short-range fermion- 
fermion interaction, broadly assumed to be a contact po- 
tential depending only on the s-wave scattering length a 
of the form 



V{v) 



Sir)- 



(1) 



This model potential, known as Fermi Gontact Inter- 
action (FGI), has been used for calculating, for example, 
the weak interaction in nucleon systems and for de- 
scribing the Darwin correction in the theory of hydrogen- 
like atoms 

FGI has been used successfully to describe Fermi sys- 
tems (cold atomic gases and nuclear matter in particu- 
lar) in the mean-field approximation, that is within the 
framework of the Bogoliubov-de Gennes (BdG) equa- 
tions [llllHZ]- However, the self-consistent solution of the 
BdG gap equation using FGI with an external potential 
diverges, leading to an infinite pairing [5]. This diver- 
gence can be renormalized as shown in Refs. [5VllO|. even 
in the experimentally and theoretically relevant unitary 
limit (a — > oo) |11H13) . where the behavior of a fermionic 
system displays universal properties. However, in this 
latter case, the Hartree-Fock (HF) term appearing in 
the BdG equations diverges as well, but this divergence 
is usually neglected [U |S] . 

Recently, few research groups have investigated the ef- 
fect of the finite-range of the potential on the proper- 
ties of dilute Fermi systems, therefore going beyond the 
FGI description. In these approaches the authors intro- 
duce an effective interatomic potential characterized by 



two length scales which allow to fix independently the 
scattering length and the potential range [TIHTB] . Previ- 
ous efforts to calculate the energy spectrum of non ideal 
Fermi p3] and Bose |I8] gases with short range potentials 
have been limited to small densities by using perturba- 
tive expansion to the second order in the parameter kpfo, 
where /o is the scattering amplitude, within Green's func- 
tion formalism. The importance of ultracold ®Li gas in 
the unitary regime goes beyond the milestone experiment 
on superfluid pairing by Ketterle |19| as it represents a 
prototype of other strongly interacting fermions, such as 
neutron matter [20 and dense quark matter j21j . 

In this paper, we show how to include density effects 
in a mean-field description of dilute Fermi gases. We 
achieve this goal by using the on-shell T matrix of the 
pair-potential as the appropriate description of the in- 
teractions and of their effective range. We will show 
that this approach removes the unphysical divergences 
obtained with a J-like potential, and allows us to investi- 
gate - in a mean field approximation - density effects in 
various observable quantities characterizing dilute Fermi 
gases. As a case-study, we solve the BdG equations [4, 22] 
for a homogeneous system of ^Li atoms in the unitary 
regime, focusing on the definition of the effective radius 
of the multichannel scattering potential, obtained from 
ab-initio calculations. At variance with approaches based 
on the momentum expansion of the s-wave phase shift, 
our theory predicts the same value of the potential range 
for both of the Feshbach resonances present in this sys- 
tem. 

This paper is organized as follows: in Sec. |Tl] we de- 
rive the principal equations. We present our approach 
where the contact potential is seen as a limit of sepa- 
rable potentials, and the BdG equations are written in 
terms of the T-matrix. This formalism is applied to ®Li 
in Sec. HI where we also describe how we obtained the 
interatomic potential using ab-initio calculations. The 
results regarding finite-range effects in this system are 
reported in Sec. ITV) 
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II. THE T-MATRIX AS THE BASIC BUILDING 
BLOCK OF BDG 

A. The contact potential as a limit of separable 
potentials 

In general terms, one can write a two-body interaction 
potential as 

V = J dridr2dr[dr'^ iPHri)^\r2) x 

F(ri,r2;r;,r^)V'(r'2)V(r'i) (2) 

where '0^ (r) is the creation operator of a particle at po- 
sition r. If, the interaction is translationally invariant, 
then the matrix elements in the previous equation fac- 
torize as: 

F(ri,r2,r;,r^) = w(ri -r2;r; - r^) x 

2 2 )■ ^-^^ 

Furthermore, a local potential is characterized by the 
condition 



i;(r; r') = (f){r)S{r — r'). 



(4) 



where (j){r) is the diagonal part of the potential and we 



have defined r = ri — r2 and r' = r'j 



We call a 



potential separable if there exists a function /(r) so that 



;(r;r') = /(r)/(r') 



(5) 



From the above definition sit is clear that the contact 
potential is both local and separable. In the usual renor- 
malization schemes, FCI is obtained as the n — >■ oo limit 
of a a series of local potentials 0n(r)- In this paper we 
take the complementary route and define the contact po- 
tential as a limit of separable potentials /„(r). In our 
approach the contact potential is seen as the limit of 
a series of projectors onto a uni-dimensional manifold. 
Therefore we write it as a limit of separable potentials 
with progressively small radius: 



F(r,r')= lim /„(r)F„/:(r') 



(6) 



where is a term of this sequence (to be specified below) 
and /n(r) is defined in term of a normalized arbitrary 
function /(r) such that /n(r) = Vri? f{nv). Our idea is to 
find the sequence Vn so that the T-matrix corresponding 
to the potential in Eq. ^ describes the proper scattering 
length a of the system. In doing so we will exploit the 
properties of the on-shell T-matrix at zero energy. Its 
fundamental role in our approach is further made clear 
by identifying the radius of the effective range of the in- 
teractions using the expansion of the on-shell T-matrix 
at zero momentum. This definition of the effective range 
differs from the usual one, which is based on the value of 
the first order correction to the s-wave phase shift 



k cot(5o = ^ xflofc^ + 

a 2 



(7) 



In our approach we expand the on-shell T-matrix at zero 
energy, T+(0), around k = obtaining: 



T,(0;k,k')--'™ -'^^ 



(8) 



For k, k' w 0, c is related to the range of the on-shell 
T+ (0) matrix in the k-space or, by Fourier-transform, to 
the inverse of the range of the T+(0) in the configuration 
space. Therefore, we define the effective range of the 
potential as: 



To 



cm 
TT\a\' 



(9) 



In order to specify the sequence Vn appearing in Eq. ^ 
we remind that the scattering length is related to the on- 
shell T-matrix via the equation: 



Ana 
m 



lim(fc|T+|fc) 
fc-i-O ' ^' ' 



(10) 



and that the Lippmann-Schwinger (LS) equation T 
V + VGo{E)T can be written as 



TiE) 



1 



V-^-Go{E)' 



(11) 



where Go{E) is the free-particle Green's function at en- 
ergy E. Expanding Go{E) in the momentum space one 
obtains: 



(12) 



,2 

where K is the kinetic energy, E = and it can be 
shown that Qi — —m\{k — 0|/)p/(47r). In terms of the 
elements of the sequence /„ defined above, one has: 



(/n 



E-K + ie 



\fn 



1 



k 



iQi 



1 fc2 

2^V 



(13) 



and therefore Eq. (11) becomes 
(k|T„+|k') = 

(k/n|/)(/|k7n)/n3 



Qq 



(14) 



By choosing 



Vn = 



Qo 



Qi 

an 



(15) 



one obtains the correct scattering length in the limit n 
oo since, for k, k' w 0: 



lim (kiT 



|k') 



1 (k/n|/)(/|kVn) 

l/Vn-Qo/n^ 



4:7: a 



(16) 
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With the same choice, the hmiting value of the matrix 
elements of the potential are: 

{■k\V\k') = hm (k|/„)F„(/n|k') = 



lim 

n— j-cxD 



(k|/n 



(/n|k') 



lim 



1 (k/n|/) (/|k'/n) 



71— J-oo 77, "-^ 



Qo 



Qi 

an 



= 0. 



(17) 



Therefore, within our approach r+(0) is a (5 function for 
regular functions and has a different behavior for not 
regular ones. At variance with Huang pseudo-potential, 
for any g{r) ^ 1/r at the origin one obtains, by using the 
definition in Eq. (15): 



lim 



lim (k|K|5) 

n—>-oo 

1 (k/n|/) n^f\g) 



(0|/)(/|ff) 



(18) 



while, as by Eq. (17 1, for regular functions the potential 



matrix elements are zero. 

Notice that, from Eqs. ([6]), (10), and (11) we have, 
denoting by |0) the zero- momentum state and setting 

9n = (/n|Go|/„), 



47ra 



- lim (0|/„)- 



1 



9n 



-(/n|0). (19) 



In the n — > oo limit, the sequence of Vn^{x\fn) is 
proportional to the S function, and the matrix element 
n^gn tends to infinity. As a consequence, one has 



lim 



Vn 



0. 



(20) 



B. The Bogoliubov— de Gennes equations for a 
separable potential 

The Hamiltonian of a system of fermions in an external 
potential U (x) and interacting with a two-body potential 
V is, in the second quantized formalism, 

H = Y^ {T^j + Uij) a\aj + ^ ^ Vu^n, a\a\a„,aj, (21) 

ij il,jm 

where a| and are the creation and annihilation op- 
erators for a complete set of single particle states. The 
indices of these operators describe all the relevant quan- 
tum numbers. 

A mean-field solution of the ground state of the Hamil- 
tonian (21) can be obtained by introducing an effective 
Hamiltonian Hcs, where the two-body potential V in 



Eq. (21 ) is substituted by an effective one-body potential 



ij 



where the as-yet unspecified matrices Wij and A^- are 
determined by requiring that the average values of H 
and H^ff — fiN are as close as possible to each other. In 
the evaluation of the average value of H, one uses the 
Hartree-Fock-Gorkov (HFG) factorization of the two- 
body density matrix which, in the case of fermions, is |24j 

{alajarnaj) = {alaj) {ala„i)- {ala„i) {alaj) + {alal) {aja„i) . 

(23) 

The self-consistent equations for Wij and turn out 
to be: 



A, 



l,jm ^;,mj) {O'l^m ) (24) 

(25) 



and one is left with a two-body effective Hamiltonian, 
which can be put in the form 



H, 



off 



En 



(26) 



where Eq is the ground state energy. The sum over k is 
restricted to those states where > 0. The new set of 
fermionic operators bk and bj, are given by the Bogoliubov 
transform 

at = {^Mjbj + v*jb^j^ 

3 



a. 



(27) 



and satisfy the relation 



exp 



(28) 



1 



although in the following we will be concerned, for the 
sake of conciseness, with the T — >■ limit. 

The equations determining the coefficients Uij and 



appearing in equation (27) are the well known 



Bogoliubov-de Gennes equations 4 : 

[(Tjfe + Uik + - tJ'5ik) Ukj + AikVkj] 
[{Tik + Uik + Wik + tJ.5ik)* Vkj + A*j,Mfcj] 



k 



eiUi 



k 



(29) 

which have to be solved self-consistently with the defini- 
tions of Wij and A^- given in equations (24) and (25 1, 
respectively. 

In our approach, where the FCI is seen as a limit of 



separable potentials, the HF term of Eq. (24 1 turns out 



to be zero. In order to see that, we remind that the Latin 
indices in the previous equations were a short-hand nota- 
tion to indicate both spin and spatial degrees of freedom. 
For example, one has 

(a|a„i) = Qai3{v,r') = E / '^ar/r, x)v^,,(r', x) (30) 



4 



where £>a/3(r, r') is the one-body density matrix, with the 
Greek letters indicating spin degrees of freedom. If we 
further assume that the pair potential is given by 



then the HF term is given by: 

T4^„fj(r,r') = lim ^(5(r-r')y / dx 

7 

War,(r,x)'u^^(r',x)] , 



r + r' R + R' 

^ 2 

(31) 



(32) 



from which, using Eq. ( pO| , we obtain that the HF term 
is identically zero. Notice that the result of this deriva- 
tion does not depend on whether the Fermi gas is free or 
confined. This result also shows that neglecting the HF 
term in the BdG equations, as usually done in the the- 
oretical mean-field treatment of Fermi gases, is indeed a 
consistent choice. 

Moreover, these results allow one to cancel out diver- 
gences in the self-consistent equations, e.g. in the pair- 
ing function, without ad-hoc renormalization procedures, 
and clarify the short-range nature of the potential in the 
sense of the limit to the 5 function, and extend the cal- 
culations to finite densities. 

To test the ability of this finite-radius approach to solve 
the divergence problem, we apply our scheme to the so- 
lution of the BdG equations for an homogeneous system 
of strongly interacting fermions at unitarity. The solu- 
tion for ultracold atoms at unitarity using FCI requires 
a nontrivial renormalization procedure j8] . 



Equation ( 33 1 can be rewritten using the Lippmann- 



Schwinger equation (11), so that the pairing function is 



determined by the normalized interaction (embodied in 
the T matrix) instead of the "bare" interaction described 
by the potential V^. 

This is achieved by rewriting the LS equation as 
V = {1 — VGn{E))T and considering the quantity 
A - VGo{E)A. Using Eq. we have the equalities 



A - VGqA 



= -VQ - VGoA (34) 
= ~V{Q + GoA) (35) 
= -{l-VGo)T{Q + GoA) (36) 



from which, assuming that (l — VGo{E)) is invertible, we 
get 



A = -r(Q + Go(£;)A) 



(37) 



which is an equation for the pairing function involving 
the T matrix instead of the "bare" potential V. The so- 



lution of equation ( 37 ) has of course to be determined self 
consistently with the solution of the BdG equations ( 29 1 . 



Introducing the indices a and /?, denoting internal de- 
grees of freedom of the atoms, Eq. (37 1 reads: 



A 



ap,p 



T+;^f},^'f}' {E; p, p') 



V 1 



E — Ea' — Epr — 



In the previous expression the value of the parameter 
E can be chosen arbitrarily and we have used this free- 
dom to have a better convergence of the self-consistent 
solution. 



C. Scattering view of the Bogoliubov— de Gennes 
equations 

In our approach, we rewrite the BdG equations at zero 
temperature in terms of the on-shell T-matrix and then 
perform our limit procedure, retaining Eq. ^ as a defi- 
nition of the effective range of the interaction potential. 

The equation for the pairing function A defined in 
Eq. ( 25 ) is readily seen to be 



A„^(r,r') = -y" dKdR' ^ 14/3,a'/3'(r, r'; R, R') 

a' f3' 

Y J dx Mc<'7(R-,x)t;^,^(R',x) 



-VQ 



(33) 



where in the last equality we have formally written the 
double integral as a "matrix product". In the case of 
the FCI potential, the equation as it stands plagued by 
a ultraviolet divergence. 



III. AB-INITIO CALCULATION OF THE 
MULTICHANNEL T-MATRIX 

We use Eq . ([38| to study an ultracold gas of ^Li atoms 
at unitarity |25), where the scattering length is infinite, 
and the observables, notably the pairing gap and the 
chemical potential, are proportional to the Fermi momen- 
tum of the free gas T. Such a regime can be obtained by 
exploiting the Feshbach resonance of a ^Li spin-mixture. 
Experimental measurements [2^ \T7\ show two Feshbach 
resonances for a magnetic strength of 543.28 G (narrow 
resonance) and 822 • • • 834 G (broad resonance) , respec- 
tively. 

The atom-atom collisions are described by an Hamil- 
tonian of the form 



H = T + H,i + Vm + Vb, 



(39) 



where T is the kinetic energy of the nuclei, Hd is the 
electron Hamiltonian (including electronic kinetic energy. 
Coulomb potential and spin-orbit coupling), and Vhf is 
the hyperfine interaction. We have also included a term 
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FIG. 1. PES for the singlet (dots) and triplet (triangles) 
states in ''Li-^Li scattering. 



Vb describing the interaction of the atoms with an ex- 
ternal magnetic field. 

In the usual Born-Oppenheimer approximation, the 
wave equation for the relative motion of a pair of ultra- 
cold atoms at zero angular momentum can be written 
as 

(40) 

where the label a denotes the total electron and nuclear 
spin states of the colliding atoms. Va{r) is the internu- 
clear potential obtained using the ground-state electronic 
wavef unction for the given value of a, and Vi^ap{r) de- 
scribes the coupling induced by the hyperfine interaction 
term, T4if, and the interaction with the magnetic field, 
Vb. 

From the solution of Eq. (40) in the E' — > limit, one 



can obtain the s-wave scattering length a of the system 
from the asymptotic values of the wavefunction, that is 



(41) 



In our theoretical analysis, the ^Li-^Li pair potentials 



appearing in Eq. ( 40 ) as a function of the external mag- 
netic field have been calculated by using Configuration 
Interaction with single and double excitations from the 
reference Hartree-Fock ground state, taking into account 
fine and hyperfine structure terms [2H]. Slater determi- 
nants have been built by expanding the molecular or- 
bitals in a cc-pVQZ basis set of atomic-centered Gaus- 
sians. The electronic singlet and triplet potential energy 
surfaces are reported in Fig. [l] using dots and triangles, 
respectively. 

In Fig. [2] we plot the calculated s-wave scattering 
length vs magnetic field for the lowest energy hyper- 
spin doublet open channel. The positions of the narrow 
(543.25 G) and broad (834 G) resonances are in good 
agreement with the experimental data j27] . 
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FIG. 2. S-wave scattering length vs applied magnetic field 
for the hyperspin states (/, m/ ) = (1/2,-1/2) and (1/2, 1/2) 
in ''Li-^Li scattering. 



The T-matrix elements appearing in Eq. ( 38 1 have 
been computed by means of multichannel scattering the- 
ory |28j . To perform this calculation we have included 
relativistic terms (of the order of 2 x 10^^ a.u.), which 
mix the singlet and triplet states, and we have multiplied 
each of the ab-initio curves of the potential by a single 
parameter fixed to reproduce the values of the Feshbach 
resonances. 

Since the open and closed channels are coupled via 
the hyperfine interaction, the problem is very complex 
if one takes into account the full on-shell multichannel 
T-matrix in the gap equation (38 1. Therefore, while the 



calculation of 7+(0) has been performed in the multi- 
channel space, we projected the T-matrix on the lowest 
hyperspin-doublet open channel to solve the BdG gap 
equation. This assumption is less drastic than using the 
interaction potential projected onto the open channel. 
In fact, the projected multichannel T+(0) matrix retains 
short range interactions with the closed channels and this 
is the reason why we prefer to use T4.(0) rather than V 
in the self-consistent equation ( 38 ) . 

The on-shell T_|_(0) is represented over 10^ equally 
spaced grid points, while a grid of 200 x 200 /c-points 
according to a Gauss-Chebyshev quadrature over the in- 
terval [0,1] has been used in the momentum space, to 
obtain a convergence below 10~^ a.u. At unitarity, we 
find that the values of the effective range ao, which is 
usually defined in term of the s-wave phase shift (see 
Eq. ([7|) are —600 a.u. and 80 a.u for the narrow and 
broad resonances, respectively. The dependence of oq on 
the resonance position and its large negative value in the 
narrow resonance rule out oq as a measure of the effective 
range of the screened potential at unitarity. 

On the other hand, by using Eq. (|9| as a definition of 
the effective radius, we obtain for rg values of 26.5806 a.u. 
and 26.5756 a.u. for the narrow and broad resonance, re- 
spectively, with a difference of less than 0.02 %. We con- 
sider this result as a further indication that the definition 
of the effective range based on the small-momentum ex- 
pansion of the T-matrix is reasonable and gives consistent 
results. 

Our methodology provides directly the values of the 



6 



extended to finite temperatures using Eq. (28). As 
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FIG. 3. Dimensionless gap function (a), chemical potential 
(b), and critical temperature (c) at unitarity for the narrow 
(dots), broad (triangles) resonances and for the 5-like inter- 
action (dashed lines). 



relevant quantities characterizing an ultra-cold Fermi gas 
- such as the gap, the chemical potential or the transition 
temperature - as a function of the effective range rg. 
Their behavior is discussed in the following section. 



IV. FINITE-RANGE EFFECTS IN THE 
UNITARY LIMIT FOR ''LI 

The ratios ^/Ep and A/Ep at unitarity as a function 
of kpro are plotted in the boxes (a) and (b) of Fig. 3j 
respectively. We observe that the behavior is exactly the 
same for the two resonances, as expected from the fact 
that we obtain the same value of the effective range for 
both the narrow and wide resonances. 

These quantities are increasingly deviating from the 
curve of the (5-like interaction (dashed horizontal line) as 
the density is increased. 

Using Quantum Monte Carlo (QMC) calculations, 
Gezerlis et al. [^Hl [^5H5T] obtained universal values at 
unitarity of A/Ef « 0.475 ± 0.075, fi/Ef « 0.4 ± 0.075, 
while Chang et al. f32} and Bulgac et al. [33] obtained 
A/Ef « 0.6 ± 0.02 and A/Ej « 0.58 ± 0.02, fi/Ej « 
0.44 ± 0.01 and n/Ej « 0.38 ± 0.06 respectively It is 
not surprising that our calculations yield results that are 
consistently higher than those provided by QMC, since 
the BdG equations are based on a mean-field approxima- 
tion of the fermion-fermion interaction. The BdG results 
are therefore very close to the BCS limit of A/Ef w 0.68 
and ^J./Ef « 0.6. 

Although the derivation of the main equations of our 
approach has been carried out in the T — > limit for the 
sake of conciseness, our theory can be straightforwardly 



consequence, we can calculate the transition temperature 
Tc - defined as A(Tc) = - as a function of the parameter 
krQ. This quantity, which is a very important parameter 
for any superfluid system, is reported in box (c) of Fig. [s] 
The ratio Tc/Tf has been previously estimated 0.494 at 
unitarity for a homogeneous system [MJ [3S] . Note that 
our theory predicts the same transition temperature at 
unitarity for both resonances, and it deviates up to 30 % 
of the corresponding FCI value by increasing the density 
of the system. 

Our results show a significant density dependence for 
A, fx and Tc, in contrast with the predictions obtained 
using FCI for which a constant value is obtained for all 
these quantities. We estimate that the largest density 
that we have investigated (corresponding to kpr^ = 0.6) 
is about one order of magnitude larger than that achiev- 
able in modern harmonic traps ^19j. Nevertheless, this 
value is likely reachable by experimental apparatus in 
a nearest future. At this high density, one might ex- 
pect that three-body recombination could lead signifi- 
cant atom losses. However, we estimate that at equilib- 
rium and for near threshold conditions 3-body recombi- 
nation will play a little role even at such increased density 
regime, since the kinetic energy gain in the recombina- 
tion process (10"^"^ a.u.) is very small compared to the 
Fermi energy (8 x 10~^ a.u.) [55] , 



V. CONCLUSIONS 

To summarize, we have proposed a new definition of 
contact potential based on the the radius of the on- 
shell T-matrix at zero energy. This analysis underpins 
a short range, rather than a (5-like, model of the screened 
fermion-fermion interaction and rules out the scattering 
length and the usual effective radius as relevant param- 
eters to describe dilute quantum gases at unitarity. Fur- 
thermore, we have shown the ability of our approach to 
cancel out naturally the divergences arising with the use 
of the contact potential. Finally, the application of this 
new theoretical approach to the self-consistent solution 
of the BdG equations for ultracold ^Li homogeneous gas 
in the unitary regime has been discussed. 
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